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Abstract: We discuss finite regions of the deconfining phase of a confining gauge 
theory (plasma balls/kinks) as solitons of the large N, long wavelength, effective La- 
grangian of the thermal gauge theory expressed in terms of suitable order parameters. 
We consider a class of confining gauge theories whose effective Lagrangian turns out 
to be a generic 1 dim. unitary matrix model. The dynamics of this matrix model 
can be studied by an exact mapping to a non-relativistic many fermion problem on a 
circle. We present an approximate solution to the equations of motion which corre- 
sponds to the motion (in Euclidean time) of the Fermi surface interpolating between 
the phase where the fermions are uniformly distributed on the circle (confinement 
phase) and the phase where the fermion distribution has a gap on the circle (decon- 
finement phase). We later self-consistently verify that the approximation is a good 
one. We discuss some properties and implications of the solution including the sur- 
face tension which turns out to be positive. As a by product of our investigation we 
point out the problem of obtaining time dependent solutions in the collective field 
theory formalism due to generic shock formation. 
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1. Introduction 

The study of black holes using the AdS/CFT was initiated by Witten's observations 
that black hole spacetimes in AdS$ are holographic duals of the deconfinement phase 
of the four dimensional SU(N) gauge theory on the boundary of AdS^ [1]. Various 
subsequent studies explored this connection to discuss the dynamics and phases of 
the gauge theory and the physics of black holes [2] - [26]. 

In a recent paper [15] two of us (PB, SRW) discussed the small Schwarzschild 
blackhole-string transition in AdS§ x S 5 by relating it in the gauge theory to the large 
N [29,30] transition, and its smoothening by non-perturbative effects. An important 
feature of studying the Schwarzschild blackhole in AdS^ x S 5 is that the blackhole 
is uniformly spread on the S 3 at the boundary where the gauge theory lives. This is 
consistent with the fact that the scale invariance of the dual gauge theory is broken 
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only by thermal boundary conditions. Consequently one can discuss the blackhole- 
string transition in terms of a unitary matrix model, where the unitary matrix (the 
thermal order parameter) is uniformly spread over S* 3 . 

In this work we would like to focus on confining gauge theories which have two 
length scales: the confinement scale A and the temperature T = f3~ l . The relevant 
order parameters in such theories have a spatial variation on the scale of A -1 . The 
gravity duals of these theories have blackhole solutions which are localized on the 
boundary. It has been argued in [14,25], that their holographic dual corresponds, 
in the large N limit, to a localised region of the deconfinement phase. This object 
has been called the plasma ball in [14], and it has a mass and a lifetime of o(N 2 ). A 
qualitative gauge theory discussion in [14] uses a balancing of positive surface tension 
and negative pressure inside the plasma ball, to argue for its existence. 

There is no doubt that it is important to study the plasma ball and its dynamics. 
Besides its utility for the physics of gauge theories at finite temperature, it is one 
more concrete laboratory for testing and studying various conundrums presented 
by blackholes [27]. The fact that the blackhole dual is localized on the boundary 
provides a greater handle on studying the horizon and what lies behind it. 

Before we begin to make headway into an understanding of these problems, we 
need to have a dynamical handle on the plasma ball in the gauge theory. This is 
a standard hard strong coupling problem. Here by strong coupling we mean large 
t' Hooft coupling A = N g\ M . One natural strategy is to use numerical techniques. 
However a direct numerical approach is also difficult without developing a formalism 
within which we can ask the right questions. 

We present a partial answer to this question in this work. We will discuss the 
plasma ball as a large N soliton which can be discussed in terms of various order 
parameters which distinguish between the confinement or deconfinement phases of 
the gauge theory. In order to do a concrete calculation we will focus on a concrete 
model that was discussed in [14] in which an interpolating solution was found between 
two bulk solutions of type IIB string theory: the AdS soliton [4] and a blackbrane. 
Both solutions are asymptotically R 2 x Si x Sg, where Si is the thermal circle of 
radius (3 and Sg is a Scherk-Schwarz spatial circle of radius 2tt. The corresponding 
gauge theory is a Scherk-Schwarz compactification of Af = 4 SU (N) gauge theory, 
on R 2 x Si x Sg . The relevant and natural order parameters of this gauge theory are 
the holonomies of the gauge field around Si x Sg. In fact for technical reasons we 
will compactify R 2 to a Scherk-Schwarz cylinder, so that the Euclidean spacetime of 
the gauge theory is R x S T x Sq x S\. The radius of is chosen larger than that 
of the Si and S} \ 

We discuss the effective action of the gauge theory in the long wavelength ex- 
pansion defined by the confinement scale A. The effective action, in the axial gauge 

1 This additional compactification of the boundary does not disturb the bulk solution 
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along the non-compact direction x, is a one dimensional model of three unitary ma- 
trices U(x), V(x) and W(x) corresponding to the zero modes of the Wilson loops on 
x Sg x S^. Using the fact that we are working with a confining gauge theory of 
adjoint fields which are all short ranged (of the order of A" 1 ) one can integrate out 
V(x) and W(x) to arrive at an effective action involving the single unitary matrix 
U(x), which has the general form 



where A -1 is the confinement scale, and f(U) and g(U) are gauge invariant functions 
of U. f(U) and g(U) contain the information that the gauge theory has a first order 
confinement /deconfinement phase transition. 

It is possible to discuss soliton solutions of the general multi-trace model using 
the Hamiltonian formulation together with the method of dealing with multi-trace 
operators developed in [15]. However in order to exhibit a solution we simplify the 
effective action even further and present the soliton (plasma kink) solution. It turns 
out to be the motion of the Fermi surface of the many fermion problem that is equiv- 
alent to the matrix model in the SU(N) invariant sector. This solution interpolates 
between the confinement and deconfinement phases and has energy density peaked 
at the phase boundary. 

In our investigations we realized that it is imperative to use the 2 + 1 dimensional 
phase space formulation of the classical Fermi fluid theory. The collective field for- 
malism, which is a hydrodynamical description in 1 + 1 dimensions inevitably leads 
to shock formation and singularities. It is not clear whether a finite energy density 
soliton solution can be obtained within collective field theory. The shocks are spuri- 
ous singularities due to the collective field description which correspond to the folds 
on the Fermi surface, which are inevitable. 

The plan of the paper is as follows. In section || we describe the two bulk 
geometries- the AdS soliton [4] and the black brane solution, for which an interpolat- 
ing domain wall solution was constructed in [14]. In section 0, we present a qualitative 
discussion as to how one can arrive at an effective description of the thermal gauge 
theory in terms of the holonomy matrices around the various cycles of the boundary, 
starting from a four dimensional gauge theory compactified on Scherk-Schwarz cir- 
cles. For technical reasons we will be working with a gauge theory compactified on 
two Scherk-Schwarz circles. One can have two dual effective descriptions, in terms 
of either the Polyakov line or the Wilson loop over the spatial cycle. We present the 
general class of such effective matrix models. In the following sections we will be 
working with a particular matrix model belonging to this class. This model can be 
discussed in terms of an exact fermionic description [28,30-35]. We shall also discuss 
the collective field equations [36] and indicate that their solution develops shocks in 
finite time. 
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In section ^ we discuss the phase structure of the model. This model has two 
stable phases: the confined and deconfined phases, and it undergoes a first order 
phase transition at a particular temperature. Later in the section we construct 
the soliton (kink) solution which interpolates between the two phases at the phase 
transition temperature. We then discuss some of the properties of the solution, in 
particular the surface tension of the soliton is discussed. We also present the localised 
(in one dimension) soliton solution at temperatures greater than the phase transition 
temperature, which approaches the confined phase in the two ends, and discuss some 
of it's properties. 

In appendix |A| we show that starting from the confined phase of the theory, 
where the density of eigenvalues of the Polyakov line is uniform, we reach the clumped 
eigenvalue distribution only asymptotically, and never in any finite time. In appendix 
[B]we discuss the relation of the shocks formed in the collective field theory description 
to the formation of folds in the Fermi description. 



2. Plasma ball in the large N gauge theory and dual black 
holes 

A plasma ball is a localized spherically symmetric bubble of the deconfining phase of 
a confining gauge theory. In [14] using the AdS/CFT correspondence, their existence 
was inferred by exhibiting a bulk solution that interpolates between the AdS Soliton 
[4] and the black-brane solution. The AdS soliton (AdSS), is given by the metric, 

ds 2 = L 2 a'(e +2u (dr 2 + T 27T d6 2 + duJ 2 ) + —L-du 2 ) (2.1) 

T 2 ir{U) 

where, 

T 2 .( M ) = l-(i(d+l)eT (d+1) (2-2) 

In this paper we will be working with d = 3. The coordinate 9 is periodic with peri- 
odicity 27r, and r is the angular coordinate along the thermal circle of the Euclidean 
theory, with periodicity r — ► r + j3, and the u>i are the two non-compact coordinates, 
while u is the radial coordinate. The boundary topology is R 2 x S]. x S], where S\ 
and Sg are the thermal and spatial cycles respectively. From the expression for Ti^, 
one sees that the spatial circle shrinks to zero size at a finite value of u. 
The black-brane (BB) geometry is given by the metric 

ds 2 = L 2 a'(e +2u (T dT 2 + d6 2 + duo 2 ) + -^—du 2 (2.3) 
with Tp(u) = 1 — (-A~(d+ l)e u )~( d+1 \ This metric continued to Lorentzian signature 



has a horizon. Notice that when (3 = 2n, the two metrics 2.1 and |2.3| are simply 



obtained from one other by interchanging the thermal circle with the spatial circle. 
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Since geometrically there is no difference between the two, the free energy of the two 
configurations must be the same at this temperature. For j3 < 2ir, the free energy of 
the BB geometry dominates the path integral while for (3 > 2n, the free energy of 
the AdSS geometry is dominant. In [14] a domain wall solution which interpolates 
between these two solutions was constructed. Clearly such a domain wall solution 
exists only for (3 = 2ir when the free energy of the two phases is equal. The domain 
wall is independent of one of the non-compact direction and in the other non-compact 
direction the BB and AdSS geometry are asymptotically reached at the two ends. 

These solutions can be incorporated within the IIB string theory by compact- 
ifying on S 5 , with the five-form RR flux turned on. This would then have a dual 
boundary description in terms of the Scherk-Schwartz compactification of the M = 4 
SU(N) SYM theory on a spatial cycle, with thermal boundary condition on both 
the cycle Si and Sg. The gauge theory lives on R 2 x Si x Sg. At f3 — 2n clearly the 
two circles are identical and can be interchanged. 

The above discussion suggests that a ball of large but finite radius of the de- 
confined plasma can occur as a solution to the finite temperature effective action of 
the gauge theory, at a temperature slightly above T c . At T = T c there exists a kink 
solution interpolating between the confined and deconfined phases. 

3. Gauge theories on R 2 x x Sj 

From the AdS/CFT correspondence, these bulk geometries- the AdSS geometry and 
the BB geometry correspond in the thermal gauge theory to the confinement and 
deconfinement phases respectively [1] . These phases are characterised by the expec- 
tation value of the Polyakov line, which is the trace of the holonomy around the 
thermal circle, 



Wi are the two non-compact coordinates and the 6 is the angular coordinate along the 
spatial circle, while V denotes path ordering. In particular, the expectation value of 
tr U vanishes in the confined phase while in the deconfined phase it takes a non-zero 
value 2 . Similarly one can define the holonomy around the spatial cycle Sg. 



Since the role of the two circles are interchanged in the two bulk geometries, it follows 
from the AdS/CFT correspondence, that trV = in the deconfined phase, and it 

2 This basically reflects the fact that a quark in the fundamental representation of SU(N) has 
infinite free energy in the confining phase and finite free energy in the deconfined phase 




(3.1) 




(3.2) 
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is non-zero in the confined phase 3 . At /3 — 2n, because the two geometries are 
identical under the interchange of the thermal and spatial circles, the effective action 
in terms of V should be identical to the one in terms of U. Later in this section we 
will qualitatively argue as to how one can arrive at an effective action in terms of 
both U and V and then in terms of either U or V, starting from the four- dimensional 
gauge theory. 

Since we will mainly be interested in the solution which interpolates between 
the confinement and deconfinement phases as a function of one of the non-compact 
direction, it should be possible to find the one dimensional kink solution in an effective 
one-dimensional unitary matrix model. In order to realize this in a gauge theory at 
large N, it turns out to be convenient to work with R x 5* x Sg x 5*, where the 
is the spatial circle, obtained by compactifying a noncompact direction previously 
labelled by the coordinate w 2 - We introduce the holonomy along the spatial cycle 



This would correspond to replacing one of the non-compact directions of the bulk 
geometry that we discussed earlier, with a circle without changing the solution. 
Henceforth we shall set the noncompact direction w\ = x. 

3.1 Effective action in terms of the Polyakov lines and Wilson loops 

The bosonic part of the action of the general gauge theory will be written in terms 
of the gauge degrees of freedom Ai, A T , A$, A a as well as the scalar fields $j which 
transform in the adjoint representation. Here A\ corresponds to the gauge field in the 
non-compact direction and we can choose the axial gauge A\ = 0. These fields are 
in general functions of (x, 9, r, a). Since the Scherk-Schwarz compactification breaks 
supersymmetry, the fermions are massive and the scalar fields get mass at one loop 
from quantum corrections. They can therefore be integrated out from the quantum 
effective action. Fourier expanding the gauge fields in all the circles and integrating 
out all the higher KK modes around every circle, we get a effective theory in terms 
of the zero modes of the fields: A®(x), Ag(x), A° a (x). 

This effective theory in terms of the zero modes is gauge invariant, and therefore 
we should be able to write it down in terms of the zero modes of the the holonomy 
matrices U,V and W. From now on we will use the notation U, V, W, to denote the 
zero modes of the above holonomy matrices. 

The effective action will be a function of all possible gauge invariant operators. 
The gauge invariant operators are constructed out of the Zn invariant products of 
the polynomials of U, V and W and their covariant derivatives, D X U, D X V, D X W, 

3 This reflects a gluon condensate in the vacuum of the gauge theory [37]. 





(3.3) 
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and are of the form Uitr(U li V mi W Pi (D x U) ni ...), where the exponents k, rrii, Pi, rij,etc 
are integers, such that the sum of all the exponents J2i h + m i +Pi + n % + ••• = 0. In 
the gauge A\ — 0, the covariant derivatives are the same as the ordinary derivatives. 
At sufficiently long wavelengths we neglect the higher derivative terms which are 
suppressed by powers of the confining scale A^ 1 . 

Depending on which of the holonomy matrices condense, there will be three 
phases in the gauge theory. In the BB phase XxU 7^ 0, trV = and XxW = 0. In 
the other two phases one of the spatial holonomy matrices V or W will get expec- 
tation values, while the expectation value for the other two vanish. However we are 
interested in an interpolating solution between the black brane and the AdS soliton, 
and not in the transitions involving all the three cycles. If we choose the radius 
R(S^) > R{Sl), R(Sq), at the temperature of interest, then from the supergravity 
solution it follows that the cycle never shrinks and corresponds to (W) = in 
the gauge theory. We can therefore put W = in the effective three matrix model 
to once again obtain a two matrix model. The action for this will in general be very 
complicated, with all terms that are allowed by gauge invariance. It will contain 
words of the type tr(U ni V n ' 2 U n:i ....), and also derivative terms. The general action 
in the long wavelength expansion will be of the form, 

S ef f = A' 1 J dxh(U, V)tr\d x U\ 2 + f 2 (U, V)tr\d x V\ 2 + (3.4) 

f 3 (U, V)tr(d x Ud x V^) + f 4 (U, V) + h.c 

where, A -1 is the confinement length scale, fa's are gauge invariant functions of 
arbitrary polynomials of U, V and (3 with appropriate factors of N. At (3 = 1/2tt, 
when the size of the two cycles are equal, the effective action will be invariant under 
U V. Integrating over either the U or the V we will get a single matrix model in 
terms of V or the U matrix. 

Since the theory is confining and has a mass gap, we can integrate out the V 
matrix, without worrying about infrared divergences, and we will be left with a 
model, given by 

S = A -1 J dxf{U)ti{d x Ud x U ] ) + g(U) (3.5) 

where again A -1 is the confinement scale, and f(U) and g(U) are temperature de- 
pendent, gauge invariant functions. As in ( |3.4j ) we have neglected all the higher 
derivative terms in the action, which are suppressed by powers of the confining scale. 
Equivalently we could integrate out U to arrive at a matrix model of V. 

In the sequel we will mainly study the soliton solution of the simplest of this 
class of models, given by 4 . 

S = A' 1 J dxNtr(\d x U\ 2 ) + £\trU\ 2 (3.6) 

4 This model has previously appeared in the discussion of 1+1 dimensional gauge theories [38] 
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Here we will assume that £ > which ensures the existence of a first order phase 
transition at some value of £. By rescaling x — > A _1 x, we can remove the explicit A 
dependence from the above action to get, 

5 = J dxNti(\d x U\ 2 ) +£|tr£/| 2 (3.7) 

where x is now given in units of A -1 . Hence forth we will be using this form of the 
action 5 . 



4. Analysis of the one dimensional matrix model 

In this section we will analyze the phase diagram of the unitary matrix model given 
by the action of the form (|3.6| ). The matrix model described by the action (|3.6| ) can 
be discussed using two methods. One is to use the collective field theory techniques 
as was done by Jevicki and Sakita [36] . This is basically a collective field description 
in 1 + 1 dimension. The Hamiltonian is written in terms of the density p(9, x) and 
velocity v(9, x) = dgll(9,x), where n is the canonical conjugate of p. The p(9,x) 
field is the eigenvalue density field constructed out of the matrix U, 

+ 0O 

p(9,x) = Pn(x)e 2 ™ ne (4.1) 

n=— oo 

where p n = -^tr ([/"■). For example, from the matrix model described by equation 
( |3.6|) , we get the following collective field Hamiltonian, 

H cf = /^(^ + ^)-e|pi| 2 (4.2) 
This Hamiltonian, gives rise to the following set of fluid dynamical equations, 



^^ + ^M*,9M*,fl)) = (4-3) 

Here 9 is a periodic variable defined in the range [— tt, tt] and a; is a variable defined in 
the range (— oo, +oo). The collective field approach is only valid for solutions which 
are spatially uniform,(for which v(x,9) = and J^p(x,9) = 0). The spatially non- 



5 Therefore all the quantities we calculate later in the text like the surface tension of the phase 
boundary of the soliton, for example, will be given in units of the confinement scale. 
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uniform solutions generically develop shocks in finite time, after which the collective 
field equations are not valid. 6 

A correct (and exact) way to analyze the model Q3.6| ) is to rewrite the model as 
a theory of interacting fermions [33] with the Hamiltonian (where the 'x' direction 
is identified with the Euclidean time). 

H= / ' d9^{9)dfy{9)-£\ I ' d9e ie ip(9)^{9)\ 2 (4.4) 



In the large N limit the fermion system will be classical and one can use the 
phase space density, U(p,9,x) such that, 



f —d9U(p,9,x) = 1 

J 27T 



(4.5) 



If a phase space cell is occupied then U(p,9,x) = 1 or else U(p,9,x) = 0. Hence 
U(p, 6, x) satisfies the relation 7 

U(p,6,x) 2 =U(p,6,x) (4.6) 

The Hamiltonian written in terms of the phase space density is, 

r/ ^ dpd9^-U{p,9,x) -£| [ dpd9e w U(p,9,x)\ 2 (4.7) 



N 2 J 2 

In terms of U(p, 9, x), the density and velocity p(9, x) and v(9, x) are, 

p(9,x) = [ $-U(p,9,x), v(9,x) = - [ p-pU(p,0,x) (4.8) 
J 2vr pj 2it 

In the appendix we will further discuss the relation between the phase space 
and collective field theory approach and we will interpret the shock formation as the 
formation of folds on the Fermi surface. Hence the shock singularities are artifacts 
of the collective field approach and are resolved by a more accurate treatment. 

In the following sections we will analyze the solutions of the fermionic Hamilto- 
nian (|4.7|) . We will start by describing the spatially uniform solution (phases of the 
theory) and then describe the non- uniform interpolating solution (plasma kink). 



As discussed in more detail in appendix 0, this phenomenon can be understood from the 



underlying fermionic theory. Infact if we change x — * ix and v — * —iv in equation (4.3), we get 
the inviscid Burgers equation with a source term. In [41], it has been shown, using the method of 
hodograph transformation, that the source free version of the Burgers equation develops shock in 
finite time. 

7 The relation (|4.6|) is true only at large N. At finite N, U satisfies the relation, U * U = 
cos jjyidgdpi — d p dgr)[U(p 7 9)U(p' ,Q')]\ p i =p fi>=e = U [35], which reduces to equation (4.6) at large 
N. 
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4.1 Spatially uniform solutions 

Here we analyze those solutions where the density of the eigenvalues of the matrix 
U is uniform over the direction x. In this case the location of the Fermi level will 
also be constant in x. Classically p can always be chosen to be an even function of 9. 
Then the potential in the equation (|4.7|) becomes, £(f dplf cos 8U(p, 9, x)) 2 . In the 
Hatree Fock approximation, the phase space evolution equation for a single particle 
is, 

9 = p (4.9) 
p = —2^pi{x) sin 9 

Where pi(x) = J dp^- cos 8U{p, 9, x ) and 9 = 4-8, p = 4-p. For a spatially uniform 
solution, pi is independent of x and we can integrate the above equations to get, 

p 2 = 2(E + 2£ Pl cos&) (4.10) 

where E is the energy of the particle. Therefore for a particle on the Fermi level, we 
have, 



p± = ±y/2(E f + 2Zp 1 cos6) (4.11) 

where p± correspond to the upper and lower branches of the Fermi level. Conse- 
quently 

V2 



p(9) = — y / E f + 2£p l cos9. (4.12) 



7T 



One has to satisfy the normalization condition given in equation ( |4.5| ) and the self 
consistency condition for p\, which effectively solves Ef in terms of £ and p\ 

J d8^^E f + 2£ Pl cos 8 = 1 (4.13) 

y/2 , 

d6 — cos 9^ Ef + 2^pi cos 9 = p x 



Depending on whether \-^^\ < 1 or \^^\ > 1, the integrals in equation (|4.13| ) will be 



evaluated between the limits [—9 ,9 ] } with 9 < tt, or over the full range [— 7r,+7r]. 
The former case corresponds to the gapped phase, as p{9) = outside [— 9 , 8 ].), 
while the latter case corresponds to the ungapped phase) [30]. 

One can study the different static phases of the model, by solving the self- 
consistency and the normalization conditions given in equation ( |4.13| ) simultaneously. 
This is hard to do analytically, but can be studied numerically. However it would be 
useful to have an understanding of the various phases as extrema of the potential in 
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terms of p. This potential can be obtained from the Hamiltonian given in equation 
Q4.71) , using equations (14.81 , |4.11|) to integrate over p. We then obtain, 



H = j d6 l -pv 2 + V{[p\) (4.14) 

where, 

V([p})= [ d6^--£\ Pl \ 2 (4.15) 



6 

The potential of the model is actually a function of the infinitely many Fourier 
modes of p. Note that the static phases are all of the form given by the equation 
( f4.12|) . It is therefore useful to parametrize p by 



P 



\ 



^2a n cos{n6) (4.16) 



n=0 



With this parametrization, the uniform phase solution is given by a® = j-, and all 
other a n = 0, while the gapped phase corresponds to a n = 0, for n > 1 and a , a x 
taking appropriate values. With this parametrization, the potential will be a function 
of the a n . Since all the phases of the theory lie in the plane given by a n> i = 0, it will 
be enough to restrict to this plane. We therefore parametrize ft± by the following 
form. 



p ± = ±V2(£ + 2£Cicos#) (4.17) 



We determine E in terms of G\ by the normalization condition ( |4.5|) . Then 
substituting this in the expression for the potential, the potential becomes a function 
of only one parameter C\. Then we can numerically calculate the potential given by 
the equation([4.15D as a function of C\ (see figure ffl). 



We now summarise the key points from our analysis of the phase structure of 
the model in consideration. 

• At low enough values of £, there is a single phase where p{6) = or p± = ±|. 
Here C\ — 0. This is the uniform phase of the eigenvalue distribution. 

• At £ = £ n = 0.227 there is nucleation of two phases for which piff) is no more 
a constant. Both the phases have a gapped eigenvalue distribution. One phase 
is unstable (//) and the other is stable {III)- 

• The first order phase transition between the phase / and phase III occurs at 
f = & = 0.237 and C x = 0.4408, E = 0.1711. 

• The phase I becomes locally unstable at £ = £2 = -25 

• At £ = £3 = 0.23125, and C\ = 0.3336, phase II has a gapped to ungapped 
transition, this is the point of the third order GWW phase transition. [29,30] 
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4.2 Spatially non-uniform solutions: plasma kinks 

In the previous section we have analyzed the phase structure of our model. In 
particular we saw that at £ = 0.237, the two stable phases (the confining and the 
deconfining phases) of the model have the same free energy. In this section we will 
first describe an interpolating domain wall type solution from the deconfined phase 
to the confining phase, at this value of £. Later in the section we will also construct 
a localised soliton solution which reaches the confined phase for large values of 

The confining phase is described by a constant Fermi level which is given by the 
following equations in phase space, 

P± = ±\ (4-18) 
While in the deconfining phase, the Fermi levels were given by, 

p± = ±^2(£ + 2fCicos0) (4.19) 



Therefore we are looking for solutions in which the Fermi level evolves from ( 4.1S|) 



to ( 4.1SQ . In terms of the geometry of the Fermi level it is a evolution from a band 



like to an ellipsoidal structure fig. 0. 

In terms of p, the solution has the property, 



p(0,x)->— , x-^-oo (4.20) 
lit 



V2 



p(9,x) — > — y E + 2£pi cos6*, x — > oo 

7T 
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Figure 2: A schematic picture of the Fermi levels. 

Now in general the Fermi level will be described by the vanishing of some implicit 
function f(9,p,x) = 0. In the static case, 

f(p, 9) = (p + - y/E + 2^C lC os9) (p_ + y/E + 2Zd cos 9) = (4.21) 

In the general case f(p,9,x) is not of this simple form and may have more roots. 
This corresponds to the case where the upper and lower Fermi levels develop folds 
and become multi-valued in 9. 8 

As each point in the phase space satisfies the equation, 9 = p,p = V'(8), one can 
derive the time evolution of the function / to be, 

d x f+pd e f + V'(9)d p f = (4.22) 

It would be interesting to try and solve the above equations numerically as a bound- 
ary value problem. We have not been able to do this. Instead we take a variational 
approach to the problem, and make a simple but reasonably accurate ansatz for the 
Fermi level. We will now summarise the main steps of the analysis. 

• We choose an ansatz for the Fermi level similar to the form in the static case, 
f(p, 9, x) = (p+ - p(0, x) + v(8, a:))(p_ + p(0, x) + v(8, x)) = (4.23) 
with p given by, 

/2 

p = — y /E(x)+2^C 1 (x)cos6 (4.24) 

7T 

where the E(x), C\{x) are functions of x. This would be a good approximation 
if the E(x), C\(x) are slowly varying functions of x. What we are doing in 
effect is to approximate the actual solution by a two Fermi surface solution 
throughout the evolution of the system, always given by the two curves p = p±. 
Therefore p± are of the form, 

p± = ±V2 V E i x ) + 2 £CiO) cos 9 + v(9, x) (4.25) 

8 In fact as is shown in appendix the folds are inevitably formed no matter what Fermi level 
configuration one starts with. 
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d 

dx 



E(x) is determined in terms of C\(x) by the condition (4J5) or equivalently 

J dO^^E(x) + 2£d{x) cos 9 = 1 
We determine v(9,x) by the continuity equation, 

[u{p,9,x)dp— = (4.26) 

J 27T 

The solution of the continuity equation is given by, 



Next, substituting this form of p{9,x) and v(9, x) back into the Hamiltonian 
and performing the 9 integral, we get, 

H = C[ 2 K(d) - V{Ci) (4.28) 

where C[ = 4-Ci(x). Hence the whole problem is reduced to a quantum 
mechanical problem of Ci(x). The function K(Ci) and V(C\) are determined 
numerically, and K(Ci) is positive and non-zero. Along the propagation in x 
the quantity TC is conserved. This conservation law is used to determine the 
relation, 



The above equation is integrated numerically to obtain Ci( function of 

x. Knowing C\(x) enables us to determine the phase space density U(p,8,x). 
The plot of C\ ( function of x is shown in figure (tfl). It should be noted 

that the soliton rises slowly but approaches the other end relatively fast. This 
follows from the asymmetric nature of the potential. 

It is important to check for the self consistency of this ansatz. This can be 
done by substituting the pi(x) obtained from our ansatz into the single particle 
equations and see how they evolve in x under this p\. One can then compute 
the p\{x) obtained from this exact evolution at each instance of x, which we 
denote by Pi(x) and compare with pi(x) obtained from the ansatz. If p\{x) 
were an exact solution, then one would get Pi(x) = pi(x). This is checked 
numerically. We started with 50 x 50 particles uniformly distributed over the 
phase space region p G [— |, ^],9 G [0, 2-71"]. This gives us the band like Fermi 
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Figure 3: Plot of C\(x) (green), p(x) (red) and free energy density (blue, not in scale) 

level in figure |2|. We study the evolution of the individual particles under the 
driving force 2£px(x) and calculate the Px(x) from the phase space distribution 
of the particles. We present the plots comparing the two values of pi(x), in 
figure £|. 

• One may also look at the snapshots of the phase space particles. In figures |9] 
and 10 we have presented two snapshots taken at x ~ 11.9 and at x ~ 11.6. We 
see from the plots that the system is driven to the gapped phase configuration 
to a good accuracy. The phase space snapshot at the later value of x matches 
very well with the expected Fermi distribution in phase III at £ = 0.237. This 
means that we are indeed reaching very near to the phase III. 

We also find that during the evolution of the Fermi sea, folds are formed on 
the Fermi level. As we discuss in the appendix |B| this is inevitable. However 
the area under the folds is a small fraction of the area of the full Fermi surface. 
This shows that our ansatz of a Fermi level with no folds, is self-consistent. 

One also sees from the phase space plots that, as discussed in appendix [A], 
p(0, x) 7^ for all x. 

• If we continue to plot the evolution of the phase space particles for long times, 
we will see that the value of Pi(x) will start falling from it's value in the gapped 
phase, and the particles will disperse away from the ellipsoid as the system will 
move away from the gapped phase. This happens because even though the p\ 
we obtain from our ansatz drives the system very near to the gapped phase 
starting from the uniform phase (as is evident from the phase space plots), it 
does not take it exactly to the gapped phase, since no matter how good the 
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Figure 4: Plot of p\{x) (green) and p\{x) (red) with x 

ansatz is it is not the exact solution 9 . If we continue to the evolve the system 
this error will start accumulating and the system will again disperse away from 
the gapped phase. This problem would not occur if we could do the exact 
numerical simulation for the soliton in the phase space as a boundary value 
problem with value of p\(x) fixed at both ends. 

An important quantity that we can determine from our solution is the surface 
tension. The surface tension in general could either be positive or negative at the 
phase boundary. However, for lagrangians with positive kinetic terms, which is true 
in our case, the surface tension also turns out to be positive. 

In one dimension surface tension is defined as the total free energy of the soliton, 
which in turn is the total action for the soliton. Hence the surface tension a is, 
(see [39]) 

/+oo 
dx(y{C x {x))-V vacuum ) (4.30) 
-oo 

This quantity at £ = £i = 0.237 is numerically calculated to be, a = 0.0027. 
4.3 Localized soliton- plasma ball 

In the previous section we constructed an interpolating kink solution for £ = £i. 
For £ between £i and £2, the two minima corresponding to phase I and phase III 
have different free energies (figure [T) and in particular, the minima corresponding to 

9 This is clear since in the correct solution folds are always formed no matter how small. 
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Figure 5: Plot of V(C\) showing the bounce solution below. 



phase J(Ci = 0) is a false vacuum. In this case there exists a soliton solution which 
is localized in the x direction, and which goes to G\ = at both x — > ±oo [40]. 

Such a solution has a simple interpretation in terms of a particle in real time 
moving in a potential —V(x). From the conservation of the Hamiltonian ( |4.28| ), it is 
obvious that if we start from C% = at x = 0, the solution never reaches phase 
It will bounces from a finite value of Cb and comes back to the phase / again, where 
Cb is determined by the relation V{Cb) = V(0). In Fig H we present a schematic plot 
of — V(Ci) and the bounce solution. 

As before, one can construct such a solution numerically (see figure ^|). This 
solution has a natural interpretation as a bubble of deconfined plasma within the 
confined phase. The plots shows two interesting trends. The first one is that the 
width and height of the soliton both increases as £ — > £i = 0.237 from above. The 
second one is that as £ — > £2, the height of the soliton decreases, but the width of 
the soliton also increases. Hence width of the soliton comes to a minimum at some 
value of £ between £1 and £2. 

One can define the width w of the localized soliton as a measure of the spread 
in x over which the value of G\ drops to a specified fraction C* of it's maximum C&. 
As £ — > £1 the localized soliton becomes the semi- infinite soliton discussed in the 
previous section and consequently the width of the soliton goes to infinity. It would 
be interesting to calculate the change in the width of the soliton with £ as £ — > £1. 
In this limit, Cb almost reaches Cni- The equation of motion is given by, 

2K(C 1 )C'I + 2K'{C 1 )C[ 2 = V'id) (4.31) 
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Figure 6: Plot of C x {x) as a function of x at £ = 0.245 (violet), £ = 0.24 (red), £ = 0.238 
(green) and £ = 0.2375 (blue). 



where C[ = £d, K'(d) = ^K(d), and similarly V'{d) = ^V(d)- Expand- 
ing K{C\),V{C X ) around C\ = Cm, and using the fact that V'iCm) = 0, and C[ 
will be small and negligible near C\ = Cm (because Cm is a turning point), we get 
from equation( |4.31| ) 

-^SC, = AiCuJSd (4.32) 



where Sd = Cm - d and A = ^jg^j . 

Using the boundary conditions, 5Ci(0) = (C/// — and ^5Ci(0) = 0, one can 
solve the above equation to obtain, 

5d = (C IH - a)(cosh(v / Ix) (4.33) 

If we define B = Cb — C#, then the width w; is given by, 

1 + — = cosh(^/Aw) (4.34) 

Cm — Cb 

Since Cm — Cb — > as £ — > £i, it follows that in this limit the leading £ dependence 
of Cxrj — Cb will be of the form Cm — Cb ~ (£ — £i) a , where a could be any real 
positive number. Putting this dependence back into the above equation, and solving 
in the w — > oo limit, we get, 

woc-log(£-£i) (4.35) 
Hence we see that the width of the soliton diverges logarithmically with £ — £i- 

5. Conclusion 



In this paper we have a presented a o(N 2 ) soliton solution of a confining gauge theory 
which interpolates between the confining and deconfinement phases separated by a 



first order phase transition. The soliton is a solution of the large N, long wavelength 
effective action of the gauge theory expressed in terms of the thermal order parameter 
(Polyakov line). The general three dimensional effective Lagrangian would have to 
contain higher derivative terms to support a soliton solution and this would make the 
problem technically very difficult. However, in the present work we have analyzed a 
simpler one dimensional example. We have presented a qualitative discussion on the 
possible connection of this model with a higher dimensional confining gauge theory 
which has a gravity dual. The soliton that we have found numerically is a finite 
region of the deconfinement phase (plasma kink/ball) with a positive surface tension 
at the phase boundary. The free energy density is also a smooth function every where 
in space. 

Even though the soliton solution is obtained in a thermal gauge theory formu- 
lated in Euclidean spacetime it is reasonable to expect it to be a static solution in 
Lorentzian spacetime at finite temperature. 10 This fact can be inferred by observing 
that the bulk solution can be analytically continued from Euclidean to Lorentzian 
spacetime. Given these facts it is tempting to identify the phase boundary as dual to 
the horizon of the blackhole. A more precise understanding of this correspondence 
will enable us to explore the structure of blackholes, especially 'inside the horizon' 
and address very directly the persistent question of the blackhole singularity. 
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A. Analysis of the clumping in the eigenvalue distribution in 
finite time 

In this appendix we will prove that if we give a small perturbation around phase 
/, p{0,x) never becomes near the point 9 = 0, at any finite x. Let us solve the 
equations of motion for individual phase space points near 9 = 0. Near 9 = we can 
make the approximation, sin# ~ 9. The equations of motion can be written as, 




(A.l) 



In this case the holonomy matrix V(x) may be a more appropriate order parameter. 
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where, 



M W = (° >) (A.2) 



Here we start by approximating p x (x) with with a step function such that 

pi(x)=pi, x>0 (A.3) 
= 0, x < (A.4) 

The solution of the equation is given by the condition, 

If we look at the Fermi level given by, p±(0) = ±p , then at "time" x the position of 
the Fermi level will be, 

MX) = coshW% Pl ) (A ' 6) 

As |pi(x)| < 1 , p±(x) does not reach at any finite time. Similar result seems 
to be true for a time dependent p±. Consequently, eigenvalue density function p(6) = 
p+(0) —P-{0) is always non-zero at the point 9 = 0. Hence any gap in the eigen value 
distribution can not open in finite time. However, the solution may asymptotically 
reach a gapped phase. 



B. Shock formation in the collective field equations and folds 
on the Fermi surface 

In this section we will show that the collective field equations develop shocks in 
finite time which can be understood from the underlying phase space picture as the 
formation of folds on the Fermi surface. The collective field equations may be derived 
from a classical theory of fermions. Consider first the theory of free fermions. We 
are looking at the phase space description of this theory. The motion of individual 
phase space points are described by the equations, 

9 = p,p = (B.l) 

From the above equation we can determine the equation of motion for a particle 
on the Fermi surface to be, 

d x p + pd e p = (B.2) 
where p denotes the value of p at any point on the Fermi surface. 
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Figure 8: Fermi level 



Now if the profile of the Fermi surface is such that for each value of 6, there are 
exactly two points lying on the Fermi surface, one on the upper and lower Fermi 
level each (like in figure ([^)), then we have 

d x p± + p±dop± = (B.3) 

where p± characterize the points on the upper and lower Fermi levels respectively. 
The source free version of the collective equations in ( |4.3| ) are simply linear combi- 
nation of the above two equations (see [33]), governing the dynamics of p + + and 
p + —p-, which are proportional to v and p respectively from ( |4.8| ). 

This identification with the collective field equations is perfectly fine for a fluc- 
tuation of the form shown in the figure (0). However because of the equation of the 
motion, points of the curve which are higher, have greater velocity than the lower 
points, hence even if we start with a simple profile like that given in figure (^), the 
profile changes due to the unequal velocity of the various points lying on the Fermi 
level to a profile of the form given in figure (|j). In figure^), where the profile be- 
comes multi-valued, the identification is not as before, since there are more than two 
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values of p corresponding to the same value of x. For instance, if at a point the Fermi 
profile has a multi valuedness of the " order four" , that is there are four values of p 
corresponding to the same value of 9, then the equation for p becomes 

2irp(6) = dp U(p, 9)+ dp U(p, 6) (B.4) 

JpZ J pi 

and similarly for the equation for pv. One can easily see that one cannot derive the 
simple collective field equations in this case. Hence the collective field equations do 
not describe the dynamics of the Fermi surface at all times. 

However we can still look at the the topmost value of p as p + and the lowest 
value of p as p- . In that case the equations governing the dynamics of p + + p_ and 
p + — ]3_ are the same collective field equation throughout, but then we see clearly 
from figure (El), that the values of these variables jumps at 9 = 9 , and hence the 
9 derivative blows up at this point. This jump will correspond to the shock of the 
collective field equations. Note that the description in terms of the fermion phase 
space is always perfectly smooth since it is after all the theory of free fermions. 

In our case we are dealing with a 1+1 dimensional interacting Euclidean fermionic 
theory given by a Lagrangian of one fermionic field ^(9) 

C = j ' dO&dJf + \d e y\ 2 + 2£ J d9d9'¥{9)^{9)cos{9 - 9')¥{9')^{9') (B.5) 



These equations give rise to the equation of the form (|4.3[). The phase space argu- 
ments discussed here will continue to hold even in this case again leading to shock 
formation in finite time (see fig. 10). But the theory viewed as a theory of fermions 
will still be valid. 
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igure 9: Phase space particles (red) at x ~ 11.9 showing the match with Fermi surface 
;reen) in phase III. 



1 




-1 ' < < < < < < 1 

-5 -4 -3-2-1 1 2 



igure 10: Phase space particles at x ~ 11.6 showing shocks at around 9 ~ —0.8 and 
ps 0.8. 
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